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1. THE PROBLEM 

On a square lattice (we consider here a D=2+l 
lattice for simplicity) the exact rotational symme¬ 
try is confined to rotations that are integer mul¬ 
tiples of 7r/2. If we attempt to build operators of 
spin J from basic components that are identical 
- that is, related by a lattice symmetry - then 
there is an ambiguity: 

T)7T 

exp{iJ9) = exp{iJ'0) W9 = — (1) 

if 

J' = J + m \fN (2) 

Thus, for example, we cannot distinguish J=0 
from J=4. 

It is conventional to assign the lowest value of 
continuum J to such an operator. However even 
if this choice is correct for the lightest state of a 
given ‘J’ - see the discussion in |^ - the ambi¬ 
guity cannot be avoided once we start calculating 
several excitations of the same ‘J’. 

Of course this is a self-inflicted problem. As 
a —> 0, we know that we recover full rotational 
invariance on physical length scales. Therefore we 
should be able to construct operators that become 
as close as we like to spin J, for any J, as a ^ 
0. This will require using linear combinations of 
operators that are approximate rotations of each 
other, where the approximation is such that it 
becomes exact as a ^ 0. 

In this poster we are going to explore how well 
the simplest embodiment of this idea works. We 
shall do so by attempting to calculate the J=0 
and J=4 glueball masses in the 2-1-1 dimensional 
SU(2) lattice gauge theory. 


2. FIRST STEPS ... 

We start with some closed loop on the lattice 
that is symmetric about the x-axis. Call it 4>ax- 
There is a corresponding loop about the y-axis, 
4>Ay, which is identical in the sense that e.g. 

{(j)Ay{t)(j)Ay{^)) = {(l}Ax{t)(l)Ax{0)) Vt (3) 

If we sum these two, 4>a = <l>Ax+4>Ay, then we ob¬ 
tain an operator (jjA that we would normally call 
‘ J=0’. (As usual we always take zero-momentum 
sums of such operators.) 

We now construct a loop that is symmetric 
around the diagonal which is at 7r/4 to the x- 
axis. We choose a loop that ‘looks’ as though it 
is roughly a rotation of ipAx- There will be an 
identical loop that is rotated by 7r/2. Summing 
these two loops gives us the diagonal operator pD- 
Clearly our trial J=0 operator is going to be 

(pj^o = pA + CpD (4) 

and our trial J=4 operator is going to be 

pJ=A = pA — CpD (5) 

(since exp(iJ0)=-l for J=4 and 9 = 7r/4.) 

If pA and pz) were indeed identical (up to a 
rotation) then we would take c=l and then our 
‘J=0’ state would contain no J=4 and vice-versa. 
However they are not identical and so we will 
choose c so as to enhance whatever approximate 
identity they possess e.g. 

{pDit)pDiO)) {pAit)pA{0)) vt (6) 

We now give the simplest example of such op¬ 
erators. Let Uj be the matrix on the link in di¬ 
rection j. I suppress the site index. Then our 
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example consists of the following path ordered 
products 

c^A = Tr{U,U,UyUMul} 

( 7 ) 

and 


= Tr{U,UyU,UyUlululul} 

(8) 


with, in each case, the addition of the operator 
that is rotated by an angle of 7r/2. So the ba¬ 
sic on-axis operator is a 1 x 2 rectangular loop, 
and the corresponding ‘diagonal’ operator is an 
ordered product of 2 plaquettes. We can obvi¬ 
ously extend this to 1 x n rectangular loops and 
to diagonal operators that are ordered products 
of n plaquettes. 

This operator is only useful if (j)D is close to 
being identical to 4>a (up to a normalisation) in 
which case the operators are close to being it/A 
rotations of each other and we can form good J=0 
and J=4 combinations. In Fig.l we compare these 
operators, using ‘twice-blocked’ link matrices ||]. 
(The normalisation c in eqns(4,5) has been chosen 
so that the vacuum-subtracted correlation func¬ 
tions at t = 0 are equal.) We observe that the 
on-axis and diagonal operators are indeed ‘close’ 
to being identical and so provide at least some 
kind of starting point for distinguishing J = 0 
from J = 4. 

3. A CALCULATION 

As a starting basis of operators we use 1x2, 
1x3, 1x4 and 2x4 loops, with the diagonal 
loops chosen in the obvious way. We construct 
these loops not just out of the elementary lattice 
links but also out of blocked links (to enhance the 
overlaps onto the lightest states). In practice we 
shall use once, twice and thrice blocked links so 
that our initial basis contains 4 x 3 = 12 opera¬ 
tors. As seen above, choosing the constant c in 
eqns(4,5) so that 

c2($,5(<)®i3(0)) = ($A(t)$A(0)) (9) 

for the vacuum-subtracted operators, with t=0 
(or indeed t=a), works quite well. We will, from 
now on, choose to use t=0, and all the opera¬ 
tors we use shall be vacuum subtracted and nor¬ 
malised so the t = 0 correlators are unity. 



Figure 1. On-axis (•) and diagonal (o) 1 x 2 loop 
vacuum-subtracted correlations. 


The exploratory calculations shown here are 
obtained on a modest run of 5000 sweeps on a 
24^ lattice at /3 = 9 in the D=2-|-l SU(2) lattice 
gauge theory. 

If our operators really were J=0 and J=4 - 
which can only occur in the continuum limit - 
then they would be mutually orthogonal: 

($j=o(0)$j=4(0)) =0. (10) 

We calculate these overlaps and find that the 
worst operator (in the sense of having the largest 
such overlap) is the one based on the 1x2 loop. 
Hence we discard it from now on, so that our basis 
is reduced to 9 operators. For these operators we 
find that the overlaps in eqn(lO) are very small 
compared to the overlaps (4>j=o‘f’j=o)- This con¬ 
firms that our procedure is working reasonably 
well. 

How do we extract a mass? We know from 
positivity that 

(4>j(t)$j(0)) < (11) 

where mj is the lightest mass (or energy). 

So given a set of operators <i)j; i=l,..,N, a vari¬ 
ational estimate of the mass is provided by 

mj 


( 12 ) 
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Figure 2. t = a effective masses from J = 0 (•) 
and J = 4 (o) correlators. 


where we have defined a set of effective masses by 
= (13) 

In Fig.2 we show the values of = a) 

for our nine J=0 and J=4 operators. We use 
t=a both because the effective masses become 
rapidly less accurate if we increase t; and also be¬ 
cause we know that at large enough t the lightest 
scalar mass will dominate. That is to say, a good 
‘J=4’ operator will have two quite distinct effec¬ 
tive mass plateaux: the one at lower t will give 
toj =4 while the one at larger t will give toj=o (if 
it is lighter). As a—> 0 the first plateau will ex¬ 
tend to ever larger values oft. In this preliminary 
test of the method we do not attempt to illustrate 
such behaviour. 

In our plot we see a distinct seperation between 
the J=0 and J=4 effective masses at t=a and, ap¬ 
plying the above variational criterion, we obtain 
the following mass estimates (in lattice units): 

to(J = 0) = 0.786(12) (14) 

to(J = 4) = 1.607(27) (15) 

The J = 0 mass is close to the expected value [0 . 
Is the lightest J = 4 state a glueball or a state 
composed of 2 J = 0 glueballs in an L = 4 partial 


wave? The latter is unlikely: the L = 4 suppres¬ 
sion near the branch point of the cut should lead 
to an effective mass that is considerably above the 
energy at the branch point E = 2m(J = 0). In 
addition we expect our single trace operators to 
have a small projection on states composed of two 
colour singlets. Evidence for these expectations 
comes from earlier numerical work ||^ . Although 
all this will eventually need explicit investigation, 
in this preliminary study we assume it to be so 
and that what we have calculated is the lightest 
J = 4 glueball mass. 

A better variational calculation would be to 
work with all linear combinations of the N(=9) 
$j. This is discussed in 0. 

Of course, simple as all this is in principle, in 
practice it rapidly becomes something of a Heath 
Robinson construction if one wants to extend it 
to arbitrary J. What one wants is a sytematic, 
general procedure that does not require tailoring 
for the individual ‘bulges’ of each J. For example 
one can proceed as follows. As a —*■ 0 the den¬ 
sity of lattice sites within two concentric circles of 
radii r and r ', where these radii are fixed in phys¬ 
ical units, becomes arbitrarily dense. This will 
even be the case if we choose r ^ r' in physical 
units, as long as |r — r'|/a ^ oo. In that case we 
can construct triangular operators (constructed 
out of some well-chosen gluonic propagators con¬ 
necting the centre to 2 sites between r and r') 
that become rotations of each other for arbitrary 
angles of rotation as a —> 0. These can then be 
used as part of a systematic procedure for con¬ 
structing operators of arbitrary spin J as a —> 0. 
This approach is described in [Q . 
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